Why an ac magnetic field shifts the irreversibility line in type-II superconductors 
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We show that for a thin superconducting strip placed in a transverse dc magnetic field - the typical 
geometry of experiments with high-T c superconductors - the application of a weak ac magnetic field 
perpendicular to the dc field generates a dc voltage in the strip. This voltage leads to the decay of the 
critical currents circulating in the strip, and eventually the equilibrium state of the superconductor 
is established. This relaxation is not due to thermally activated flux creep but to the "walking" 
motion of vortices in the two-dimensional critical state of the strip with in-plane ac field. Our 
theory explains the shaking effect that was used for detecting phase transitions of the vortex lattice 
in superconductors with pinning. Some recent experiments on this subject are discussed. 

PACS numbers: 74.60. Ge, 74.76. Bz, 74.76. Db 



Experimental investigation of the equilibrium proper- 
ties of type-II superconductors can be performed only 
in the reversible region of the magnetic field (H) ver- 
sus temperature (T) plane. But often strong flux-line 
pinning prevents the measurement of equilibrium prop- 
erties. In this context, Willemin et al. JlJ recently made 
an interesting observation. Their experiments revealed 
that application of an additional small oscillating mag- 
netic field perpendicular to the main dc field leads to a 
fast decay of the currents circulating in the critical state 
of various high-T c superconductors. This effect dramat- 
ically extends the observable reversible domain in the 
H-T plane. The relaxation of the irreversible magne- 
tization in these experiments was exponential in time, 
and thus was obviously different from thermally activated 
flux creep, which leads to a logarithmic time law. Us- 
ing this vortex-shaking process, the melting transition 
of the vortex lattice in YBa2Cu30 7 _,5 crystals was de- 
tected at temperatures very close to the critical tem- 
perature T c , where the melting could not be investigated 
before. With the same shaking method, it was discov- 
ered H that the order-disorder transition in the vortex 
lattice of Bi 2 Sr 2 CaCu208 is of the first order at low tem- 
peratures, where vortex pinning usually masks the corre- 
sponding jump in the equilibrium magnetization. Thus, 
the vortex-shaking process opens new possibilities in ex- 
perimental investigation of the H-T phase diagrams of 
superconductors. Nevertheless, the nature of this impor- 
tant effect so far remained unclear. 

In this paper we give a quantitative explanation of this 
vortex-shaking effect for thin strips. It is shown that the 
relaxation is caused by the generation of a dc electric 
field by the ac magnetic field. We obtain this result in 
the framework of a quasistatic approach based on the 
standard critical state concept. Therefore, this effect is 
quite general and should occur in all type-II supercon- 
ductors, not only in high-T c materials. In contrast to the 
generation of an electric field in the case 0] when the ac 



and the dc magnetic fields are parallel, the above elec- 
tric field appears in superconductors without a transport 
current. We also point out a new way of measuring the 
critical current density j c in superconductors. 

The dc electric field in general depends on the orienta- 
tion of the ac magnetic field with respect to the currents 
circulating in the superconductor. In the following we 
shall analyze the simplest situation: A thin supercon- 
ducting strip fills the space |x| < w, \y\ < oo, \z\ < d/2 
with d <C w; the constant and homogeneous external 
magnetic field H a is directed along z, while the ac mag- 
netic field hcostut is applied along x, i.e., perpendicular 
to H a and to the currents in the sample (Fig. 1). We also 
make here the usual Bean assumption that the critical 
current density j c does not depend on the local induc- 
tion B. The field H a is assumed to be sufficiently large 
to exceed both the field of full penetration for the strip, 
|5) H p = (j c d/7r) ln(2eiy/d), and the lower critical field 
H c i, and so we may put B = hqH. 

A qualitative explanation of the shaking effect is as fol- 
lows: The currents flowing in the critical state of the strip 
generate a nonuniform distribution of the magnetic in- 
duction B z (x). The ac field periodically tilts the vortices 
in this state. However, at each point x with a nonzero 
sheet current J(x) (the current density integrated over 
the thickness d) , the tilt is not symmetric relative to the 
central plane of the strip z = 0, and during each cycle 
of the ac field, the asymmetry leads to a shift of vortices 
towards the center x = of the strip. This shift tends 
to equilibrate B z (x), and it also generates a dc electric 
field which decreases J(x). When J(x) reaches zero, the 
asymmetry disappears, and the process stops. 

It will be shown below that the generated dc electric 
field is proportional to the thickness of the sample, d. 
Thus, to describe the effect, the strip cannot be consid- 
ered as infinitely thin, and the appropriate critical state 
problem is two dimensional. However, according to the 
approach of Ref. 0, the smallness of the parameter d/w 
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enables us to split the problem into two simpler prob- 
lems: A one-dimensional problem across the thickness of 
the sample, and a problem for the infinitely thin strip. 
Namely, we first interpret a small section of the strip 
around an arbitrary point x (see Fig. 1) as an "infinite'' 
slab of thickness d placed in a perpendicular dc magnetic 
field B z (x) and in a parallel ac field h cos Lot and carry- 
ing a sheet current J(x). The resulting dc electric field 
E y = E y (J,B z , h) for the slab we then use as the local 
electric field E y {x) for an infinitely thin strip, to calcu- 
late the temporal evolution of the sheet current J{x) and 
induction B z (x) in this strip by the method of Ref. |7| . 

We begin with the analysis of the dc electric field E y 
generated by the ac magnetic field, h cos tot, in an infinite 
slab with sheet current J and with a constant and homo- 
geneous magnetic field B z . At sufficiently small to |§ , the 
problem may be considered quasistatically: We solve the 
critical state equation for the slab, dB x /dz = ^oj cr {z, t), 
at every moment of time t, and then find the shape of the 
flux lines at this moment from dx/B x = dz/B z and the 
condition that the points on each line where j cr changes 
its sign cannot move. The critical current density j cr 
is always equal to j c or —j c , and its distribution over z 
is just the well-known distribution in a slab in a paral- 
lel magnetic field and with current J. In particular, at 
t = and t = 27r/w one has j cr = j c for z_ < z < d/2, 
and j cr = —j c for —d/2 < z < z_ (see Fig. 1), while 
at t = tt/lo one has j cr — —jc f°r z + < z < d/2, and 
jcr = jc for —d/2 < z < z + where z + = — z_ = J/2j c - 
An essential finding of this analysis is that during one 
half cycle (0 < t < n/u>) every flux line turns around a 
fixed "swivel point" with the coordinate z — z + , while 
during the second half cycle (w/ui < t < 2w/u}) the line 
turns around another fixed point with z — z_. As a re- 
sult, the line "walks" along x, see Fig. 1. The shift of the 
vortices during one full cycle is 
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where h p (J) = (J c — \ J\)/2 is the field of full penetration 
of parallel flux into a slab with current J, and J c — j c d 
is the critical sheet current. Note that in deriving Ax 
we have assumed h > h p (J). Else, the x component of 
the magnetic field does not penetrate completely into the 
slab, a section of the flux line does not sway and hence 
does not move at all, and thus Aa; = 0. The dc electric 
field generated by the above shift of the vortex lines is 
E y = (lu/2tt)AxB z , and hence we arrive at 
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for h < h p (J) , 
[ft - h p (J)] for h > ft p (J) . (2) 



The described mechanism for the generation of a dc 
electric field was proposed many years ago [|| , and Eq. (2) 
is equivalent to the formula (2.8) of Ref. Equation (2) 



also coincides with the corresponding result for the slab 
obtained in our paper Q , though in Ref. 0] the dc mag- 
netic field lies in the plane of the slab. This coincidence 
H is not surprising, since in both cases E y may be inter- 
preted as generated by the transfer of the x component 
of the flux across the slab thickness d. 

In the derivation of E y we have assumed that J and 
B z do not change during one cycle. This approximation 
is justified when the relaxation time of the profiles J(x) 
and B z {x) considerably exceeds the period 2-k/lo of the 
ac field. As will be shown below, this condition is indeed 
fulfilled for thin samples, d <C w, and thus the approxi- 
mation is self-consistent. 

We now consider the temporal evolution of the pro- 
files J{x) and B z (x) in the infinitely thin strip. JtJ This 
evolution is caused by the slow "walking" of vortices to- 
ward the center of the strip and is given by the Maxwell 
equation: 
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with E y specified by formula (2). Since E y — v x B z , 
where v x is the average velocity of the vortices, the above 
equation simply expresses the conservation law of their 
number inside the strip. On the other hand, the Biot- 
Savart law yields a relation between the sheet current J 
and B z : 
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Equations (2)- (4) together with the initial condition 

J(x)\ t =o = -J c x/\x\ , (5) 

are sufficient to determine the temporal evolution of 
J, E yi and B z in the strip, i.e., the functions J(x,t), 
E y {x,t), and B z {x,t). 

To proceed further, it is convenient to rearrange 
Eq. (3), (4) as follows: We invert JTo) Eq. (4), differ- 
entiate both sides of the result with respect to t, insert 
Eq. (3), and arrive at an equation for J(x,t): 
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with E y (J) given by formula (2). The right hand side of 
Eq. (6) is proportional to (d/w)u>; this becomes evident if 
one introduces the dimensionless length x/w. Thus, the 
decrease of J / J c during one cycle is determined by the 
small parameter d/w, and the above assumption that J 
and B z are constant over one cycle is well justified. 

The numerical method of solving Eq. (6) is well elab- 
orated, @] and we use it to analyze the evolution of the 
sheet current J(x, t). In Fig. 2 we show the time depen- 
dence of the magnetic moment per unit length of the 
strip, M{t) = J_ J(x, t)x dx, at various amplitudes h 
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of the the ac magnetic field. If h/J c > 0.5, then |M(t)| 
decreases almost exponentially with time, and eventually 
the equilibrium state [|l2| is established since J — > ev- 
erywhere in the strip; when h/J c < 0.5, the equilibrium 
state is not reached, Fig. 3. At large times t, the satura- 
tion value of the current in the sample, Joo, can be found 
from the condition Ey^J^) = 0, yielding 

Joo = Jc - 2h . (7) 

In the special case h = h* = 0.5 J c , Eq. (2) gives E y oc J 2 , 
and it follows from Eq. (6) that M oc J oc £ -1 . Thus, 
we conclude the following: When in an experiment the 
amplitude h is fixed while j c increases due to a decrease 
of temperature T, then at some T one has J c = 2h, and 
at lower T the vortex-shaking will not lead to the equilib- 
rium state. In other words, the greater is j c d, the greater 
is the ac amplitude h required for the equilibrium state 
to be reached. This feature of the vortex-shaking process 
was indeed observed in the experiments. |||| 

We consider now more closely the practically impor- 
tant case h > J c /2. When J decays, then, according to 
Eq. (2), the problem formally reduces to the Ohmic strip 
problem analyzed in Ref. Q , but with the resistivity now 
associated neither with free flux flow nor with thermally 
assisted flux flow, namely, E/j = p = noi/d 2 (2h/ J c — 1) = 
const, v = lo/2tt. In this case the sheet current after some 
transient time tends to the following solution of Eq. (6): 

J(x, t) = -CJ C f{x/w) exp(-t/r) , (8) 

where C is a constant depending on the magnetic history, 
f(v) is the normalized eigenfunction with lowest eigen- 
value A of the integral equation, 

ir z J _ l u — v \1 — v z / du 
and the relaxation time r is given by 

T - 1=A ^(2h-JA 

W V Jc J 

Numerical analysis |tJ yields A = 0.6386. The odd 
function f(v) = — /(— v) normalized to J Q vf(v)dv = | 
is shown in Fig. 4. The magnetic moment per unit 
length of the strip, M = J_ w xJdx, is then M = 
— Cw 2 J c exp(— i/r), while in the fully penetrated criti- 
cal state one has M = Mq = —w 2 J c and in the Meissner 
state @ M = — ttw 2 U a . 

Thus, we obtain that at h > J c /2 the decay of the 
current becomes exponential in time, in agreement with 
the experimental data |l]] , and the spatial profile of the 
current tends to the universal function f(x/w), which 
is independent of the superconducting parameters and 
dimensions of the strip. It also follows from formula 
(10) that the rate of the decay increases with u and 



h. However, for the above quasistatic analysis to be 
valid, u> should not be too large ||, while the condition 
(tuj/2tt) > 1 leads to (w/d)(J c /h) > 1, or h < j c w. 
This puts an upper bound on h up to which the above 
formulas are applicable. 

Interestingly, with Eqs. (7) and (10), measurements of 
the amplitude h* at which = 0, or of r for h > h*, 
allow one to obtain j c at short times t ~ w" 1 where j c is 
only weakly affected by thermally activated creep. 

To conclude, we give a quantitative description of the 
vortex-shaking effect in a thin superconducting strip. It 
is shown that in the framework of the critical-state the- 
ory, the decay of the critical currents to zero, or the de- 
crease of the magnetic moment of the strip to its equi- 
librium value, is due to the generation of a dc electric 
field by an ac magnetic field applied normal to the main 
dc magnetic field. This electric field originates only if 
the finite thickness of the strip is taken into account, and 
so the theoretical description of the vortex-shaking effect 
is, strictly speaking, a two-dimensional spatial problem. 
But the small value of the parameter d/w enabled us 
to simplify the problem and to solve it explicitly for ac 
amplitude h -C j c w. When the thin superconductor has 
a different shape, e.g. looks like a square platelet, there 
are regions where the ac magnetic field is parallel to the 
sheet current. In this case, the analysis of the problem 
becomes more complicated. Nevertheless, the obtained 
results explain the properties of the vortex-shaking effect 
observed in the experiments |l|-§l and also shed light on 
the origin of similar intriguing effects |p^|-|l5| discovered 
earlier for slightly different geometries. 
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FIG. 1. Strip geometry and magnetic fields (upper inset). 
A flux line "walking" from left to right through a section 
of the strip ("slab") is shown at times tio/ir — 0, 1, 2, 3. 
The swivel points of the line are marked by circles. Here 
h/J c — 0.6, J/Jc — 0.5, yielding z + — —z- = 0.25d and 
Aa; = 0.35(/io J/B z )d. Ax is the shift per cycle, Eq. (1); x 
is measured from an arbitrary point of the "slab". The line 
shape consists of parabolas with d 2 x/dz 2 — ±/ttoJ c / B z . The 
scheme at the right shows the current profiles across the slab 
at the extremal times 1 and 2 and at some intermediate time, 
where the two arrows indicate the penetration of the cur- 
rent-inversion front. 
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FIG. 3. Profiles of the sheet current J(x, t) in the strip at 
various times t (in units to /WO, to = -Kw/cud) starting from 
Eq. (5). Top: h/J c = 0.3; the dashed line shows J^, Eq. (7). 
Bottom: h/J c = 0.7; exponential relaxation to J = 0. 
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FIG. 2. The time dependence of the magnetic moment 
M{t) of the strip at various amplitudes h of the ac mag- 
netic field: h/J c = 0.1, 0.2, ... 0.8, 1, 2. The dashed line 
(h/Jc — 0.5) separates the complete and incomplete relax- 
ation of M. M(0) = Mo = —w 2 J c is the magnetic moment of 
the strip in the fully penetrated critical state. The time unit 
to = (irw/ujd) corresponds to w/2d S> 1 cycles. 
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FIG. 4. The universal spatial profile f(v) of the relax- 
ing sheet current J(x,t) in the strip at t > r, Eq. (8), 
and the magnetic field H z generated by this current pro- 
file. The dash-dotted line indicates — H z (0). H z is in units 
J/f = CJ c exp(-t/r) = -M(t)/w 2 . At the edges H z has a 
logarithmic infinity, which is cut off by the finite thickness d. 
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